Quantum Walks with Non-Orthogonal Position States 
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Quantum walks have by now been realized in a large variety of different physical settings. In 
some of these, particularly with trapped ions, the walk is implemented in phase space, where the 
corresponding position states are not orthogonal. We develop a general description of such a quan- 
tum walk and show how to map it into a standard one with orthogonal states, thereby making 
available all the tools developed for the latter. This enables a variety of experiments, which can be 
implemented with smaller step sizes and more steps. Tuning the non-orthogonality allows for an 
easy preparation of extended states such as momentum eigenstates, which travel at a well-defined 
speed with low dispersion. We introduce a method to adjust their velocity by momentum shifts, 
which allows to experimentally probe the dispersion relation, providing a benchmarking tool for the 
quantum walk, and to investigate intriguing effects such as the analog of Bloch oscillations. 



Quantum Walks (QWs) are a widely used model sys- 
tem for transport processes. Initially introduced from 
a computer science perspective [IH5], the field has sig- 
nificantly expanded and is now largely treated from a 
physics perspective [5141 0j. In fact, "quantum walk" is 
now widely taken to be synonymous with "discrete time/ 
discrete space quantum dynamics" of a particle with in- 
ternal degrees of freedom. On one-dimensional lattices, 
a QW can always be implemented by a concatenation 
of coin operations and successive state-dependent shifts 
[TT] . Already these one-body systems are capable of sim- 
ulating various physical effects such as Anderson local- 
ization [12] or the formation of molecules [10]. In par- 
ticular, single-particle QWs are a basic building block 
in a bottom-up approach towards general-purpose multi- 
particle simulation environments 13 . Therefore, one of 
the main interests in QWs is the possibility to study key 
features of quantum dynamics in a setting which can be 
controlled experimentally with high precision. 

Experimentally, QWs have been implemented in sev- 
eral different ways, for example using nuclear magnetic 
resonance [Hj , atoms in optical lattices [TS] , trapped ions 
[SHU] or photonic systems QU [20H27] . 

In the theoretical description it is almost universal 
practice, to model the different "positions" by mutually 
orthogonal subspaces in Hilbert space. However, orthog- 
onality cannot be achieved in some proposals [281130] 
and the related experiments [16-18 , which use coher- 
ent states of a harmonic oscillator for the position states 
(Fig. [I]). In order to fit the theoretical model, it was 
necessary in the experiment to choose the step size in 
phase space sufficiently large to make these states ap- 
proximately orthogonal. 

The aim of this letter is to show that the lack of or- 
thogonality can be exploited. First of all, we give a 
complete analysis of QWs with non-orthogonal position 




FIG. 1. (Color online) nQW implemented in a harmonic- 
oscillator phase space. The positions are coherent states, il- 
lustrated by their Husimi functions, P(a) — Kola^)! 2 , for 
each \a x ) separately. The inlay illustrates their orientation 
in phase space, implementing a nQW along a line. Since the 
position states \a x ) are coherent states, they are not orthog- 
onal. The step size Act — \ai — ao\ = y/2/cr (Cf. eq. |2|) 
of the nQW determines their overlap. QWs of this type have 
been implemented with trapped ions [16H18| . 

states (nQWs) , and introduce a transformation to the or- 
thogonal case, such that all results known in that case can 
be utilized. Therefore, the overlaps between different po- 
sition states no longer need to be avoided. In experiments 
with trapped ions, one can therefore consider smaller step 
sizes and thus run the walk for more steps before the re- 
quired Lamb-Dicke approximation breaks down [17] . 

Secondly, the transformation of the nQW into a QW 
with orthogonal position states encodes the properties of 
the non-orthogonality into the initial state of the QW. 
Therefore, utilizing the overlap, one can prepare sev- 
eral interesting initial states directly, in particular those 
which are extended over several positions, such as consid- 
ered in [31] . In contrast, in the orthogonal case with a lo- 
calized initial state, the preparation requires a more elab- 
orate process, including several additional operations, 
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some of which must involve a breaking of the trans- 
lational symmetry. Such a preparation process would 
severely decrease the fidelity of the experiment. 

Furthermore, we show how the initial state can be 
shifted in momentum space by including an additional 
operation into the walk operator. This allows for the 
control of the scaling and the measurement of the dis- 
persion relation, providing a benchmarking tool for the 
QW. Finally, we use this method of momentum shifts to 
implement Bloch oscillations [32] as an example for the 
range of experiments with nQWs, which can be readily 
implemented using state-of-the-art technology. 

Throughout the paper, we relate our theory to the 
trapped- ion setting (Fig. [I]). However, it applies to arbi- 
trary unitary nQWs in any dimension. 

We consider the Hilbert space % — £2(2) <& C 2 , with 
£2(2) the position space and C 2 the coin space. The (nor- 
malized but not orthogonal) position states \ot x ), with 
x £ Z, form a basis of ^(Z). The coin states are given by 
the er z -eigenstates |c+) and |c_), where <r x ,y,z denote the 
Pauli matrices. We assume the initial state of the nQW 
to be localized at the origin, that is, po = |«o) («o I ® poo 
with poo = |c+) (c+| [33]. One step of the nQW is given 
by the application of the walk operator W = S- (1 ® C), 
which is composed of a unitary coin operator C and a 
unitary shift operator S. The latter acts as 



S\a k ) 



(1) 



In fact, starting from |ao), the subsequent application of 
S defines all other position states. Hence, their overlap 
(afc|a;) is translation invariant. Modeling the trapped- 
ion systems (Fig. [IJ, we define the overlap function 



g{k) = (a x \a x+k ) = exp(-/c 2 /cr 2 ) 



(2) 



for all x, where a determines the overlap between differ- 
ent position states. 

The probability to find the walker at position \a x ) after 
t steps is related to the projector F x — \a x ) (a x \ [T51ll7| . 
That is, 



Pt{x) 



Tr(G-p ) 



(3) 



where we introduced G = T (8> I with the Gram matrix 
r = J2 X \ a x) (&x\ for the normalization. Note that the 
normalization is independent of the step number t, be- 
cause [S, G] = 0, which can be checked using the unitarity 
of S. 

In Fig. [2] the position probability distributions of 
nQWs with two different realizations of C are illus- 
trated. The first one (Fig. [2^,) is with the coin operator 
Ce = exp (i 7r/4 a y ). This type of coin has been imple- 
mented experimentally with the initial state po [16H18) . 
We therefore refer to it as the experimental walk. The 
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FIG. 2. (Color online) Position probability distribution Pt(x) 
Q after t — 150 steps of a nQW with overlap function g(k) Q 
and (a) the experimental coin Ce and (b) the Hadamard coin 
Ch for a = {0, 1.5, 8} (green, blue, red). (For the green curve, 
only points x with P(x) 7^ are connected.) In the orthogo- 
nal case (green) the probability distributions of both types of 
walks are equal, however for large overlaps they differ signifi- 
cantly (blue, red). In case (a), the probability distribution ap- 
proaches a Gaussian shape centered at the origin of the walk, 
as a is increased. The spreading, which is still linear in the 
step number t, is vastly reduced |17| . In the Hadamard case 
(b), the probability distribution approaches a shape consist- 
ing of two Gaussian peaks centered around ±t/\/2. Thus, the 
(linear) spreading is increased, as the probabilities between 
the peaks vanish. The initial state is po = |a?o) (ao|®|c+) (c+|- 



second coin operator (Fig. |2p) is the Hadamard ma- 
trix, which can be written as Ch = ct z -Ce- The coin 
operators Ce and Ch are very similar and indeed, the 
probability distributions are equal in the orthogonal case 
(c = 0). However, in the case of large overlaps (a > 1), 
they show significantly different behaviour. 

In the following we will transform the nQW into an 
orthonormal basis and show that W generates a QW in 
that basis, but with the initial state being in a superpo- 
sition of several position states. 

The Gram matrix T gives the relation between the ba- 
sis {|a x )} and its dual basis {\a' x ) = lo^)}, fulfilling 
(a x \a' y ) = S xy [34]. This allows us to define an orthonor- 
mal basis {lea,) = T~ x /' 2 \a x )} with r -1 / 2 being hermi- 
tian. Since the shift operator S commutes with G, its 
action in the orthonormal basis is 



S\e k ) 



c± 



(4) 



That is, W also defines a QW in the orthonormal basis. 
Since the probability to find the walker in position \e x ) is 
related to the projector r -1 / 2 F x T -1 / 2 , we can transform 
Eq. ^ to 



P t (x) = Tr((\e x ){e x \®l)-W t p W 



(5) 



where the initial state amounts to 

P0 Tr(GpoG) (6) 

with po = |eo) (eo| <8> poo- In particular, po is extended 
over several position states according to the overlap func- 
tion g(k). However, the fact that po is extended does not 
already imply the properties of the nQW. As shown in 
Fig. § the experimental and the Hadamard walk show 
entirely different spreading, although they differ only by 
the coin operator. In the following we investigate the 
properties of the nQW using Fourier methods and asymp- 
totic perturbation theory [3"5W37| . 

Given a vector tfj = J2 X \ e ^l ® l^x) in ^(Z) ® C 2 , 
its Fourier transform in the conjugate momentum space 
L 2 ([-7r,7r),C 2 ) amounts to ip(p) = e* p ' x l^x)- That 
is, we consider the Fourier transformed vectors as de- 
valued functions of p. 

The walk operator W is translation invariant on li^K) 
and thus acts as a multiplication operator in momentum 
space, i.e. (Wip){p) = W(p)ip{p), with W{p) = S{p)-C 
and S(p) — exp(ipa z ). From the eigendecomposition 
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w( P ) = J2^ h(p)p M (7) 

k=l 

we obtain the dispersion relations u>k(p) an d the corre- 
sponding eigenvectors ipkip), with Pk(p) denoting the 
projector onto V'fc(p)- The eigenvectors Vfc(p) define 
Bloch waves with distinct momentum p. The role of the 
dispersion relations u>k(p) is the same as for a particle 
in a periodic potential, e.g. an electron in a solid-state 
system. It encodes the fundamental transport proper- 
ties of that system. In particular, the group velocities 
Vk(p) — &u)k{p)/&p (Fig. [3^,) determine the spreading 
behaviour of the initial state of the QW. Precisely, the 
ballistic order of the spreading (i.e. linear in t) can be 
captured by the time-asymptotic position probability dis- 
tribution Poo (9), where q G [—1,1] denotes the asymp- 
totic scaled (oc 1/t) position variable. P^ (q) can be 
computed as the inverse Fourier transform of the charac- 
teristic function 

C(A)= / dpTr(po(p).e* Ay W) (8) 

J [ — 7r,7r) V 7 

with V(p) — J2k v k{p)Pk the group- velocity operator 
|37j . The initial state amounts to po(p) = \g(p)\ 2 - Poo, 
where g(p) is the Fourier transform of the overlap func- 
tion g(k). Therefore, for each momentum p, \g(p)\ 2 deter- 
mines the influence of the corresponding group velocities 
Vk(p) to the asymptotic probability distribution. 

The group velocities of the cases Ce and Ch are the 
same, but shifted by p — tt/2 (Fig. [3^). Therefore, in the 
orthogonal case both walks lead to the same probability 
distribution, since \g(p)\ is constant in p. That is, all 
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FIG. 3. (Color online) (a): Group velocities Vk(p) (grey, dot- 
ted/dashed) of the Hadamard- walk and ||po|[(p)/(27r) of the 
initial states localized at p — (blue) and p — tt/2 (red, cor- 
responds to the experimental walk) with a — 4. (b): Asymp- 
totic position probability distribution Poo{q) for each initial 
state (blue, red). For each initial state only the group veloci- 
ties around their points of localization determine the position 
probability distribution. That is, since the blue initial state 
(a) is localized at p — 0, where Vk{0) ~ ±l/\/2, Pcc(q) con- 
sists of two peaks moving away from the origin with that 
velocity (b). In contrast, the red initial state (a) is centered 
around Vk(n/2) w 0, which leads to a localized asymptotic 
position probability distribution Paa(q). The coin part of the 
initial states is poo = |c+) (c+|. 

velocities Vk(p) occur with equal weight, leading to non- 
zero probabilities in the whole range x € [— f/v2, t/y/2\ 
(Fig. [2]). The maximal velocities Ufe(p) = ±l/y/2 play a 
special role by the formation of caustics, leading to the 
well-known peaks at x = ±t/V% [37] . 

In the non-orthogonal case, g(p) is localized at p = 
0, such that only the related group velocities occur in 
the nQW (Fig. [3^i). In the Hadamard walk, Ufc(0) = 
±l/\/2, i.e. the velocities that are also most pronounced 
in the orthogonal case, whereas in the experimental case 
Vk(0) = 0, such that the position probability remains at 
the initial position. 

Due to the small but finite width of g(p), also group 
velocities close to p = influence the nQW. Since in the 
case Ce they change strongly around p = 0, the width 
of the peak in position space increases linearly in t (See 
[T?] for numerical results). Similarly, also in the case Ch 
the widths of the two peaks in position space increase 
asymptotically linearly in t, but at a much smaller rate 
(Cf. the finite widths of the peaks in fig. [3]d). 

In the following we introduce a method to shift the dis- 
persion relation in momentum space and thus to change 
the group velocity of the nQW. For a momentum-shift 
of the amount of O, we apply after each step the oper- 
ator 1 (£> i?(8) = 1 ® exp(i8<7 z ). Due to the identity 
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R(Q)S(p) = S(p + O), this is equivalent to a nQW with 
the effective walk operator 



w e (p) = s(p + eyc. 



(9) 



The time evolution is then determined by the group ve- 
locities Vk(p + Q)- Thus, using the experimental coin Cg, 
it is possible to achieve the spreading of a Hadamard walk 
by including the operator R(—ir/2) into W. In fact, since 
Ch = o z Cei the momentum shift with O = — tt/2 com- 
pensates for the avfactor, up to a complex phase. 

The momentum-shift method allows for the experi- 
mental determination of the dispersion relations ujk(p) 
by implementing the nQW with Wq for several values of 
G [—if, 7r] and determining the scaling of the position 
probability distribution for each. This is particularly im- 
portant if the walk operator is not exactly known, e.g. 
due to experimental imperfections and in the regime of 
a high number of steps. 

In semiconductor-superlattices, the driving force for 
Bloch oscillations is implemented by a static external 
electric field, leading to a linear drift of the electrons' 
momentum. The periodic band structure causes the os- 
cillatory behaviour of the electrons, detectable by optical 
methods [35]. The analog of this motion can be imple- 
mented in QWs by applying the momentum-shift opera- 
tor R(t- AO) (modulo 2-k) at the f-th step (for every t), 
which implies a shift in momentum space by AO at each 
step of the nQW. The walker thus experiences different 
group velocities at each step of the walk (Fig. [3^,) , which 
results in an oscillating behaviour in position space, in 
contrast to linear spreading with a constant group veloc- 
ity (Fig. [I]). Note, that this effect does not require the 
non-orthogonality of the position states. However, the 
simple shape of the position probability distribution of a 
nQW (two distinct Gaussian peaks) can reduce the effort 
for detection. 

A convenient system for the experimental implemen- 
tation are trapped ions. Recently, one-dimensional QWs 
with three, resp. 23 steps have been realized in the phase 
space of the harmonic motion [151 HH] and a protocol 
for the extension to 100 steps has been proposed [FT] , 
The number of steps in these experiments was limited by 
two requirements: On the one hand, the motional am- 
plitude of the ion needed to remain small, because the 
implemented protocols [38J were designed assuming the 
Lamb-Dicke approximation |39j . On the other hand, the 
step size Aa (Fig. [IJ was chosen sufficiently large in or- 
der to minimize the overlap between neighbouring posi- 
tion states. For the implementation of Bloch oscillations, 
a small step size is favoured, leading to a significantly 
higher possible number of steps. A possible choice of pa- 
rameter values for a Lamb-Dicke parameter of 77 < 0.3, as 
in Ref. [17], is given in Fig. [4] The momentum-shift op- 
erator i?(0) can be implemented by shifting the phase of 
the driving light fields with respect to the relative phase 
of the coin states [3D]. The positions of the peaks can 
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FIG. 4. Probability density (greyscale with black: Pt(x) = 1) 
of a Hadamard-nQW (a = 14) with Bloch oscillations. After 
20 steps of the nQW (Cf. fig. [2J3), the Bloch oscillations are 
switched on with AO = n/10, such that the positions of the 
peaks oscillate with a period of 20 steps and an amplitude of 
5 positions. The Bloch oscillations are switched off after 65 
steps, a point where the group velocity is zero (Fig. [3^.), such 
that the peaks remain at their position during the remain- 
ing nQW. The expectation values of the harmonic-oscillator 
occupation-number operator N range during the oscillations 
from {Nmi n ) = 1.3 to {N max ) — 2.7, which is detectable with 
state-of-the-art trapped-ion technology [39] . 



be determined using state-of-the-art blue-sideband proto- 
cols |39j . During preparation of our manuscript we were 
made aware of a related implementation of Bloch oscil- 
lations [41] , In contrast to that approach, our method 
does not require a position-dependent coin operator and 
may therefore require less technical effort. 

In summary, the transformation of nQWs into orthog- 
onal ones allows for an intuitive understanding of their 
properties in terms of the dispersion relations. Hence, by 
the momentum-shift method it is possible to change the 
spreading behaviour, to determine the dispersion rela- 
tion, and to implement - by the correspondence to solid- 
state systems - the analog effect of Bloch oscillations. 
Therefore, the non-orthogonality can be exploited and 
need not to be avoided, leading to a higher number of 
steps and a new range of experiments with available tech- 
nology. Furthermore, nQWs can be considered for mod- 
eling transport processes in complex systems and may 
lead to a better description than previous approaches. 
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